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Thermal noise is expected to be the dominant source of noise in the most sensitive frequency band 
of second generation ground based gravitational wave detectors. Reshaping the beam to a flatter 
wider profile which probes more of the mirror surface reduces this noise. The "Mesa" beam shape 
has been proposed for this purpose and was subsequently generalized to a family of hyperboloidal 
beams with two parameters: twist angle a and beam width D. Varying a allows a continuous 
transition from the nearly-flat {a — 0) to the nearly-concentric {a — tv) Mesa beam configurations. 
We analytically prove that in the limit D —> oo hyperboloidal beams become Gaussians. The ideal 
beam choice for reducing thermal noise is the widest possible beam that satisfies the Advanced LIGO 
diffraction loss design constraint of 1 part per million (ppm) per bounce for a mirror radius of 17 
cm. In the past the diffraction loss has often been calculated using the clipping approximation that, 
in general, underestimates the diffraction loss. We develop a code using pseudo-spectral methods 
to compute the diffraction loss directly from the propagator. We find that the diffraction loss is not 
a strictly monotonic function of beam width, but has local minima that occur due to finite mirror 
effects and leads to natural choices of D. For an q = tt Mesa beam a local minimum occurs at 
D = 10.67 cm and leads to a diffraction loss of 1.4 ppm. We then compute the thermal noise for 
the entire hyperboloidal family. We find that if one requires a diffraction loss of strictly 1 ppm, 
the a = 0.9l7r hyperboloidal beam is optimal, leading to the coating thermal noise (the dominant 
source of noise for fused-silica mirrors) being lower by about 10% than for a Mesa beam while other 
types of thermal noise decrease as well. We then develop an iterative process that reconstructs the 
mirror to specifically account for finite mirror effects. This allows us to increase the D parameter 
to 11.35 cm for a nearly-concentric Mesa beam and lower the coating noise by about 30% compared 
to the original Mesa configuration. 

PACS numbers: 04.80.Nn,07.60.Ly,41.85.Ew,42.15.Eq 



I. INTRODUCTION 

The initial baseline design for the Advanced LIGO 
gravitational wave detectors [l], 0| employs Gaussian 
beams in the arm cavities. The leading noise source 
in the most sensitive frequency band of the instruments 
(~ 30 — 300 Hz) is the thermal noise in the substrate and 
reflective coating of the mirror test masses. Lowering 
thermal noise is therefore of paramount importance for 
achieving a higher event rate in LIGO. There are a num- 
ber of other detectors that are being built or upgraded 
to similar specifications. While we will choose to study 
Advanced LIGO for definiteness, our general conclusions 
should be more widely applicable to any interferometeric 
detector that needs to limit thermal noise. Some of the 
important parameters that we use are summarized in Ta- 
ble HI 

LIGO is a Fabry-Perot interferometer with four mir- 
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rored test masses. The resonant beams in the cavity mea- 
sure the position of the test masses, averaging over the 
mirrored surface, with the average weighted by the power 
distribution of the beam. Thus, the highly illuminated 
central area is weighted more than the mirror boundary 
that is left nearly dark. One way of decreasing the ther- 
mal noise is to flatten the beam so that a larger fraction 
of the mirror is in use. Motivated by this intuitive rea- 
soning, O'Shaughnessy et. al. [1, Q proposed the flat 
topped Mesa beams, which were subsequently explored 
in detail by them and others 0, S 0, S] ■ These beams 
would lower the thermal noise by a factor of approxi- 
mately 2.5 compared to the baseline design. The original 
Mesa beam supported by nearly-flat Mexican Hat mir- 
rors was found to be susceptible to a tilt instability 
This triggered the proposal of a Mesa beam supported 
by nearly-concentric mirrors (lo| . In the same paper, a 
family of hyperboloidal beams that include all Mesa and 
Gaussian beams previously considered was introduced. 
Mesa is currently the leading alternative beam design 
for Advanced LIGO, and is being studied experimentally 

In this paper we first discuss the formulation of hyper- 
boloidal beams. The "nearly-flat" Mesa is created by su- 
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perposing minimal Gaussians with generators uniformly 
distributed inside a cylinder, and the "nearly-concentric" 
by generators falling inside a cone, and passing through 
the center of the cavity. These two choices have the same 
intensity distribution on the mirrors, but the second has a 
much smaller susceptibility to tilt instability. The hyper- 
boloidal beams smoothly interpolate between these two 
cases by twisting the generators of the minimal Gaus- 
sians by an angle a. After discussing some geometric 
properties of the beams, we present a proof that Gaus- 
sian beams are a special case of the hyperboloidal beams, 
confirming a conjecture in p^ . 

We then compute the three types of mirror thermal 
noise for a variety of hyperboloidal beam s hap es, using 
a set of simple scaling laws developed in [ij, [lj| that 
simplify previous work [l^ [l^, [l3| • The first is substrate 
Brownian noise, occuring due to mechanical dissipation 
in the material; this is the least significant source of ther- 
mal noise. The substrate thermoelastic noise is caused 
by random thermal expansion. The coating also has both 
Brownian and thermoelastic noise, but these follow the 
same scaling laws so we consider them as a single source. 
The coating noise is the most severe of the three types 
for the fused-silica substrates currently planned for Ad- 
vanced LIGO. The substrate thermoelastic noise would 
dominate in a material like sapphire which has a higher 
thermal expansion coefficient. We show that the noises 
decrease with increasing width of the beam, as expected, 
and that the hyperboloidal beams have larger noise than 
the relevant Mesa beams. 

The constraint on our ability to lower the noise 
comes from the need to keep the diffraction loss small. 
Gravitational- wave interferometers must keep a large cir- 
culating power in the cavity, and so cannot allow signif- 
icant amounts of light to escape past the edge of the 
mirrors. The current design constraint used in the most 
recent papers [1, [H, [l^ is a diffraction loss of 1 part per 
million (ppm) per bounce for 17 cm fused-silica mirrors. 
The clipping approximation Eq. (jlSp indicates that the 
desired Mesa width is approximately D — 10 cm. Previ- 
ous work 0, [20| have shown that the clipping approxima- 
tion is not accurate for Gaussian beams of finite mirrors, 
and have used Fast Fourier transform simulations for ac- 
curate calculations. 

We calculate the diffraction losses accurately from 
eigenvalues of the cavity propagator using an exponen- 
tially convergent code that enables us to study the struc- 
ture of Mesa and other hyperboloidal beams as a function 
of D, mirror radius, and twist angle a in detail. We find 
that the diffraction loss is not a monotonic function of D, 
but due to finite mirror effects has anomalous local min- 
ima where the loss is significantly below what is expected 
from the clipping approximation. These minima are ob- 
served to become more shallow and eventually disappear 
as the radius of the mirror is increased. However, for 
the mirror radii and beam widths relevant for Advanced 
LIGO the finite mirror effects are important. We show 
that they can be used to increase the width of the beam. 



lowering the noise even further than previous work. 

Finally, we develop an iteration scheme to redesign the 
mirror, explicitly accounting for finite mirror effects. The 
iterated mirror is altered to match the phasefront of the 
primary eigenbeam of the finite mirror cavity, reducing 
the diffraction loss of this mode, thus allowing even larger 
beam widths to satisfy the diffraction loss constraint. 

The mathematical construction of the hyperboloidal 
beams are discussed in Sec. II, while the asymptotic 
limit of the hyperboloidal beams are derived in Sec. III. 
The thermal noise scaling laws are described in Sec. IV. 
The cavity propagator construction and eigenmode de- 
composition are presented in Sec. V, with the results, 
including finite mirror effects, discussed in Sec VI. We 
then summarize our work in Sec. VII. 



L 


3999.01 m 


Length of LIGO cavity 


Ao 


1.064 X lO"** m 


Laser wavelength 


Wo 


^\oL/2n = 2.6023 cm 


Minimal Gaussian width 


R 


17 cm 


Mirror radius 



TABLE I: Advanced LIGO Parameters 



II. CONSTRUCTION OF THE BEAMS 

The beams we study are supported by two identical 
mirrors facing each other, forming a cavity of length L. 
The mirrors are cylindrically symmetric around the opti- 
cal axis, which runs along the length of the cavity and will 
be called the z axis. The center of the cavity, equidistant 
between the mirrors, is z = 0; the mirrors are located at 
z — —zn and z — zj^, where = L/2. The transverse 
distance from the z axis will be denoted by r, and the 
angular coordinate by </>. 

The cavity is fed with laser light with wavelength Aq, 
and the distance between the mirrors is fine-tuned so that 
the cavity resonates in its fundamental mode, with a field 
amplitude U{r,z) and intensity |C/(r, In this paper 
we focus on axisymmetric modes with no (j) dependence 
in the beams. Non-axisymmetric modes are important 
for studies of the tilt and parametric instabilities [6|, |2i[ , 
but are not discussed in detail in this work. The narrow- 
est possible Gaussian mode that can exist in a cavity of 
given length is called the minimal Gaussian, which has 
the intensity distribution 

(1) 

where Wq = y/XozaJn. 

A hyperboloidal beam is the superposition of minimal 
Gaussians chosen such that the symmetry axis of the 
individual minimal Gaussians are generators of a set of 
coaxial hyperboloids. The beam family has two parame- 
ters: a, the twist angle one would have to rotate the two 
basis of a set of coaxial cylinders with respect to each 
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FIG. 1: (a) The beam intensity profile ^aS^ and (b) correc- 
tions Ha are shown at fixed D = 10 cm for different twist 
parameters a. 



other to obtain the hyperboloids and D, the radius of a 
section perpendicular to the optic axis of the outermost 
hyperboloid at the end of the cavity. In the case a = 0, 
the propagation axes are parallel and fill a cylinder of ra- 
dius D. This is the Mesa beam supported by nearly-flat 
mirrors. For a = tt the lines all cross at z = form- 
ing two cones. This configuration also generates a Mesa 
beam, but one supported by nearly-concentric mirrors. 
Varying a smoothly deforms the beam and the mirror 
shape between the two configurations. Some examples of 
the beam shapes are displayed in Fig. [T^. 

For the cavity to support the desired beams, the phase 
of the electric field of the beam should be constant on the 
mirror surface. We will focus our attention on the mirror 
on the positive z side of the cavity. The wavefront can 
be approximated by the "fiducial spheroid" , 



(2) 



For a = this is the mirror plane z = zji; for a — n, 
the fiducial spheroid is a sphere centered on 2 = 0, and 
clearly the lines which generate the hyperboloid are all 



radii of the sphere. 

There are two equivalent expressions for the field am- 
plitude evaluated on the fiducial spheroid. The first is 
the integral expression [To| 



Ua{r,Sa) = A / / d'/'o»'oexp 
Jo Jo 

^2 I J2 



■^^0 . , . 

I — TT smfflo sma 



(r^ +r;^ - 2rrocos0o) . s 
„ ^(1 — zcosa) 



2wl 



,(3) 



where A is a complex constant. 

The second is the method that we use in this paper. 
The beam is constructed as shown in [22| by an ex- 
pansion in Gausse-Laguerre eigenbeams of spherical mir- 
rors. They are closely related to the Gauss-Laguerre basis 
functions given by 



V-^CO = V2exp(-eV2)i„(C2) , 



(4) 



where Lm is the m*'' Laguerre polynomial. Then the 
Gauss-Laguerre eigenbeams are 



«'m(r',z) 



w{z) 



V2r 

w{z) 



exp 



2R{z) 



(5) 



X exp[z(/coz - (2m + l)$(z))] 



where 



w{z) = wo^/l + (z/zr)'^ , R{z) = z + Zr/z, 

$(z) = arctan(z/zii;) 

and fco — 27r/Ao. The expansion is written as 

00 



(6) 



(7) 



The expansion coefficients that result in hyperboloidal 
beams are 



cosar^p(^m + l,^ 



P(a, x) is the incomplete gamma function 

1 



r(a) Jo 
Jo 



(8) 



(9) 



Jo 



The mirror shape that supports a hyperboloidal beam 
is not exactly the fiducial spheroid. We make a correction 
h{r) so that the surface of the mirror is given by ZMir) — 
Saif) ~'h{r). The correction is chosen so that the mirror 
is located at a phasefront of the beam. We find the mirror 
surface from Ua{r,Sa) by 



h{r) = ko\Uo.ir, Sc.) - t/„(0,5„)) 



(10) 
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The shape of h{r), or the mirror itself, is generally re- 
ferred to as a "Mexican hat" , and some examples are 
displayed in Fig. [TJd. 

As expected, the beams for a and —a are identical, as 
they correspond to hyperboloids that are simply twisted 
in opposite directions. There is a duality between a and 
TT — a. This was first mentioned in [To| for all a's, then 
studied in more depth for a = and a = tt in [6| and, 
finally, analytically understood in Q. This duality ex- 
tends to several quantities. The beam intensity profiles 
are identical. The corrections to the mirror shape are 
opposite, ha{r) = —hT^-a- There are also dualities in the 
eigenvalues of the propagator [1, [l^, that we will 
not discuss in this paper. 



to evaluate the sum in p3)) . with t = icosa. The final 
result is an intensity profile 



\Uc.{r)\' 



1 



■ exp 



sma 



(15) 

The minimal Gaussian a = 7r/2 is seen to have a = wq. 
The width is symmetric under a — > tt — a, as expected 
from the duality relation 0, , and goes to infinity at 
a = or TT. This includes every Gaussian beam capable 
of resonating in a cavity of the given length. 



IV. THERMAL NOISES 



III. ASYMPTOTIC BEHAVIOR OF WIDE 
HYPERBOLOIDAL BEAMS 

It was conjectured by Bondarescu and Thorne that 
the beam becomes a Gaussian in the limit D — > oo. We 
will prove this analytically for the intensity profile of the 
beam, evaluated on the plane z — zji which would be the 
surface of a perfectly flat mirror. The intensity varies 
slowly enough with z that this will also be the intensity 
profile on the mirror to a good approximation. Our proof 
uses the expression for the beam amplitude in terms of 
a summation of Gauss-Laguerre functions. The essential 
ingredient is the realization that the expansion coeffi- 
cients take the form Am^ — [constant)™- as D oo, 
where the constant depends only on a. This allows us to 
analytically perform the summation to obtain the beam 
profile. 

In the limit a; — > cx), the incomplete gamma function 
P{a,x) = 1, giving = V2w§(- cos a)™/!?^ ^j^g 
approximation z ~ zif yields 

w{zr) ^ V2wo , R{zr) = 2zB, , ^(zr) = 7r/4 , (11) 

and the Gauss-Laguerre propagators become 

*™(f, zr) = -Lv'™(r)e-^^'"e^'^« , (12) 

where f = r/wo and the r-dependent part of the phase 
has been absorbed into (j){f). Since in the end we will 
only be interested in the intensity profile, the exact form 
of 4>{f) is unimportant. 

The expansion for Ua then becomes 

oo 2 

U^{f,ZR) = ^(-cos«)'"^(e-^)'"e-*L„(f2)e^^« 

rn— 

= (^) e-Ve'^«^(zcosa)™L„(f2). (13) 

m=0 

We now use the generating function for the Laguerre 
polynomials [11] 

f; L„.{x)t"^ = ^ exp (- Y^^) (14) 



There are a number of noise sources limiting the sen- 
sitivity of ground-based gravitational-wave interferome- 
ters. Seismic noise causes an effective cutoff in the lowest 
frequencies that can be measured. Fundamental prob- 
lems such as shot noise and radiation pressure noise, as 
well as technical issues, are important limitations on the 
sensitivity throughout the frequency band. However, the 
major contribution in the most sensitive frequencies of 
LIGO is the thermal noise in the mirrors. Reducing the 
thermal noise is the goal of this paper. 

The mirror consists of a substrate with a coating, and 
we must consider noises due to fluctuations of both. The 
substrate and coating have both thermoelastic and Brow- 
nian contributions to the noise. Thermoelastic noise is 
caused by expansions in the material caused by random 
heat flow. Brownian noise is due to the coupling of nor- 
mal modes of vibration by imperfections in the material. 
As a practical matter, both types of noises in the coat- 
ing have the same scaling law so they do not have to be 
considered seperately. In fact, for the fused silica mirrors 
now under consideration for Advanced LIGO, the coating 
noises are the dominant contribution. However, use of a 
material like sapphire, with a higher coefficient of ther- 
mal expansion, would cause the substrate thermoelastic 
noise to dominate. We will calculate all three types of 
noise in this paper. 

A set of simple scaling laws were derived by Lovelace 
[Tst in parallel with O'Shaughnessy [l3] that are applica- 
ble to beams of arbitrary shape. The noise is proportional 
to an integral depending on the intensity and an overall 
constant which is independent of the shape of the beam. 
The noises are given by 

poo 

Sn^Cn i{kfk''dk , (16) 

Jo 

I{k) = / Mkr)\U{r)\\dr , (17) 

where n specifies the type of noise under consideration, 
and I(k) is the 2D axisymmetric Fourier transform of 
the beam intensity with k the radial wavenumber. The 
substrate Brownian noise has n = 0, coating Brownian 
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and coating thermoelastic noises have n = 1, and sub- 
strate thermoelastic noise has n = 2. We are interested 
in comparing noises given by different beam shapes, so 
the overaU constants C„ are not important. The result- 
ing amphtudc sensitivity is the square root of the noise, 
and has units of meters / \/Hz . 

These scaling laws were derived for half-infinite mir- 
rors, meaning that effects of the finite radius and thick- 
ness of the mirror are ignored. For the specific mirrors 
under consideration for Advanced LIGO, the width and 
thickness of the mirror are large enough compared to 
the beam width [l3| that this should be a good approx- 
imation. The results of [l7| suggest that we can expect 
corrections of not much more than ten percent to the 
half-infinite scaling law expressions that we are using. 



V. EIGENVALUES OF THE PROPAGATOR 

The idealized picture of a locked cavity is that the mir- 
rors are perfectly aligned and a precise distance from each 
other. The beam should leave one mirror, reflect off the 
other, and when it returns it should be the same shape 
and exactly in phase. The beam will have lost some inten- 
sity due to diffraction and the finite extent of the mirror. 
In order to build up a very intense beam with a relatively 
weak laser, the beam must refiect very many times. The 
loss per half-trip (from one mirror to the other) must be 
below approximately 60 ppm. The majority of the loss 
will be due to absorption in the mirror and other factors, 
not diffraction. The commonly agreed upon budget for 
losses due to diffraction is 1 ppm. 

In previous work, the clipping approximation is often 
used to estimate the diffraction loss by calculating the 
fraction of the intensity of the beam which falls outside 
the mirror. It is given by 



Jr 



(18) 



Our numerical code is accurate and fast enough to com- 
pute the diffraction loss directly from the eigenvalues of 
the propagator, which is more accurate. 

The propagator also allows us to estimate the difficulty 
of locking the interferometer. The finesse of Advanced 
LIGO is about 1200 [2411. This sets the width of the reso- 
nance for the cavity [2^ to be about 27r/1200 ~ 0.005 
radians. If any other modes with a small diffraction 
loss have an argument within this distance of the de- 
sired mode, there will be severe problems with locking 
the cavity. For axisymmetric modes, we find that this 
is not the case and the arguments of the eigenvalues are 
well-separated enough that locking with Mesa or any of 
the other beams that we study should be no more difficult 
than locking the currently proposed Mesa. 



A. Integral Form of the Propagator 

The propagator from a single point r, cf) on one mirror 
to a point r' , (p' on the other (sec for instance _5:j) is 

/C(r,0, /,(/.') = ^(l + cos^-*^-"" , (19) 

where p is the path length between the two points, 
the angle between the cavity axis and the path, and r 
and (j) are the standard cylindrical radial coordinate and 
azimuthal angle. The cavity is very long compared to 
the radius of the mirrors so we can immediately make 
the paraxial approximation ^ = 0. The path length can 
be approximated as 

p = {(S^ir) + Sa^ir') - h{r) - h{r')f (20) 

« i + cosa(^- + -j-M^)-Mr') 

- — COs{(j) - (j)') 

since L ^ r and L 3> h{r). This expression is the path 
length used to find the phase of the propagator. In the 
prefactor of the propagator we only need the much sim- 
pler p K L m order to compute the amplitude to an 
accuracy of R^^^^^/L'^ ~ 10~^. 

We assume that the cavity is axisymmetric, which 
means that the eigenmodes of the propagator can also 
be written as eigenstates of rotation. Thus we write the 
complex beam amplitude as Vm\r.,4>) = Vm {r) e"""*^. 
To determine the modes of the cavity, we have an eigen- 
value problem given by the integral 

7^V«(r,0) = rT'dr'd<t>' /C(r,(/.,r',0')Kl^'(^',0') , 
Jo 

(21) 

where is the associated eigenvalue, and R is the 

mirror radius. Integrating over (j)' produces 

7^V«(r) = r r'dr'JCl,{ry)V^\r') , (22) 



^n{ry) ^ 



L 



k^rr' 



cxp 



ikn{-L + h{r) 



+h{r') 



cos a 



(23) 



where for convenience we have defined a "radial kernel" 
A^m- Jm{x) is the m*"^ order Bessel function of the first 
kind. In what follows, we specialize to axisymmetric 
modes (m = 0), as we are not focusing on tilt or para- 
metric instabilities which involve modes with to > 0. 
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B. Discrete Form of the Propagator 

The above integral must be converted to a discrete sum 
to be suitable for numerical computation. We choose to 
do this using a Chebyshev quadrature [H, [2^ , which is 
appropriate for the finite range of r' . We use N collo- 
cation points rfc with associated integration weights Wk- 
The integral is converted to 

7«K)«(r,) = Ef^i/CS(r„r,)rfc«;fcT/«(r,) . (24) 

This is a matrix eigenvalue problem which is easy to solve 
numerically. We will order the eigenstates by the number 
of radial nodes, the number of times that V^^*'' (r) goes to 
zero. The fundamental mode (i = 0) has zero nodes and 
as i? — > oo it limits to the hyperboloidal beam that the 
cavity is designed to support. We will show that finite 
mirror effects cause the actual eigenstate to be slightly 
different. 

We performed a convergence test where we varied the 
number of collocation points, N , and observed exponen- 
tial convergence. The diffraction loss for the lowest eigen- 
value for the a = TT Mesa beam with D = fO.67 cm, 
R = n cm changes in the 5th digit when varying N be- 
tween N = 250 and N = 500. We are using N = 1000 or 
N — 500 for all of the calculations in this paper. On a 
typical single processor laptop computer, our code takes 
~ 10 seconds for N = 500 and - 1 minute for N = 1000, 
to calculate the eigenmodes of any given cavity. 



C. Meaning of Eigenvalues 

The eigenvectors of the propagator are the field am- 
plitudes of the cavity's resonant modes. Fine-tuning the 
length L of the cavity io L~\-5L changes the eigenvalue to 
exp(— ifco(5L)7. This tuning is used to select the desired 
mode; when the argument of the eigenvalue is zero it will 
resonate in the cavity. The magnitude of the eigenvalue 
is the fraction by which the amplitude changes during a 
half-trip. We must have I7I < 1 for finite mirrors, be- 
cause some light will always be lost to diffraction. We 
define the diffraction loss in parts per million (ppm) over 
one half trip through the cavity as 



1000 



p/: = 10^(1 - 



(25) 



Advanced LIGO requires a diffraction loss per half trip of 
about one ppm. In the next section, we will discuss the 
results of our analysis, where we study the noise charac- 
teristics of the hyperboloidal beams given the constraint 
on the diffraction loss. 



VI. RESULTS 

Our goal of reducing the noises in LIGO is constrained 
by the need to keep the diffraction loss at nearly 1 ppm. 
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FIG. 2: (a) The diffraction loss for an {a = vr) Mesa con- 
figuration (a — tt) is shown as a function of D on a loga- 
rithmic plot for several different mirror radii: i? = 16 cm, 
R — 17 cm and R = IS cm. The diffraction loss computed 
numerically using Eq. (|25p (solid, dashed and dotted-dashed 
lines) exhibits local minima due to finite mirror effects. It 
can be seen that the minima get narrower as R increases and 
that they go below the values estimated using the clipping 
approximation (dotted lines), (b) The fractional difference 
1 I finite/ 1 t^l theory ~ 1 between the theoretical infinite-mirror 
beam intensity profile and the actual profile given by the 
first eigenvector of the propagator in Eq. (|22|l is plotted for 
D = 10.67 cm and R = 16, 17, 18 and 20 cm. The deviation 
decreases with R as expected. 



The hyperboloidal beams have two parameters, D and a. 
D is roughly the width of the beam, and a is the shape. 
The duality relation reduces the range of a that we need 
to consider. We focus our attention on 7r/2 < a < n 
since the intensity profiles are identical to those in the 
range from to 7r/2 but the mirrors are nearly con- 
centric as needed to decrease the tilt instability. The 
Mesa profile obtained for a = tt has the flat-top shape 
required to decrease the thermal noises. As a goes to- 
ward 7t/2 the beam becomes more rounded, losing the 
flat top and sharp decay near the edge. We have shown 
that D — > oo the beam becomes a Gaussian whose width 
is woVl + cos^ a/ sin a. For a near tt, this Gaussian is 
nearly infinitely wide, and at a = 7r/2 the beam becomes 
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FIG. 3: Diffraction ioss as a function of D is displayed for 
a = TT, a = 0.957r, a = O.QOn. It can be seen that as a 
decreases the minimum diffraction loss is lower and occurs 
for larger D. 




0.84jt 



0.887t 0.9271 

a 



0.967t 



FIG. 4: The minimum diffraction loss in ppm (represented 
by crosses) and the corresponding D (dots) are shown as a 
function of a. The solid line represents the best exponential 
fit of the form a + exp(6 + csin^ a) with a = 0.094, b = -4.34 
and c = 2.20. 



the minimal Gaussian of width wq. From the clipping 
approximation we can estimate that ai a ~ 0.2477r and 
a ~ 0.7537r, D = oo the diffraction loss is about 1 ppm. 
We do not have to consider any values of a between these 
since the widest beam consistent with the diffraction loss 
constraint would be the D = oo Gaussian. 



A. Finite Mirror Effects 

If the mirrors were infinite in extent, no light would 
propagate off of the mirror and there would be no diffrac- 
tion loss, giving eigenvalues of unit magnitude. The clip- 
ping approximation assumes that the beam is the one 
supported by infinite mirrors. This is not the case as 



a = IT 


D = 9.62 cm 


Diffraction Loss 


Phase 


1.0 ppm 


0.0 


121.7 ppm 


-1.5104 


loz.o ppm 


-U.Oooo 


334.8 ppm 


-2.6677 


7941.6 ppm 


2.2904 


45401.8 ppm 


0.8325 


a = O.OStt 


D = 9.71 cm 


Diffraction Loss 


Phase 


1.0 ppm 


0.0 


136.2 ppm 


-1.5079 


195.3 ppm 


-U.ooUZ 


302.3 ppm 


-2.6630 


7225.3 ppm 


2.2976 


44104.7 ppm 


0.8422 


a = 0.907r 


D = 11.01 cm 


Diffraction Loss 


Phase 


1.0 ppm 


0.0 


18.6 ppm 


-0.4797 


951.8 ppm 


1.3048 


3400.8 ppm 


2.7653 


5870.7 ppm 


-2.3322 


61706.3 ppm 


1.4493 



TABLE II: The diffraction loss and phase separation for eigen- 
values with losses less than 10 percent, in three different hy- 
perboloidal configurations with 1 ppm loss in the fundamental 
mode. Only axisymmetric modes are shown. 



diffraction also causes the beam profile to change. The 
propagator is a more accurate calculation because it finds 
the precise beam profile supported by the mirrors. The 
clipping approximation is typically an underestimate of 
the diffraction loss |,||. In Fig. Hi, we show that this 
is indeed usually the case. However, for some ranges of 
beam width D there is an anomalously low diffraction 
loss below the clipping approximation. 

To study this effect, we varied the mirror radius and 
computed the diffraction loss as a function of D. The 
mirrors that we study in this paper have radius R ~ 
17 cm, so we compared with R — 16 cm and R = 18 
cm. The local minimum becomes narrower and shallower 
for increasing mirror radii. This is suggestive of a finite 
mirror effect that will disappear when the mirror radius 
is significantly larger than the radius of the beam. 

Fig. [21d shows the deviation of the beam from the infi- 
nite mirror Mesa beam. A concentric Mesa {a = vr) beam 
with D = 10.67 cm is chosen. This is the location of a 
local minimum of the diffraction loss with respect to D. 
The beam intensity is computed directly from the eigen- 
vector of the propagator. The intensity is normalized by 
integrating over the mirror rather than over the entire 
mirror plane as with the infinite mirror case. However, 
this only causes a fractional error in the normalization 
on the order of 10~^. As shown in the figure, the finite 
mirror causes deviations from the infinite mirror beam; 
the plateau of the beam is less flat than expected. When 
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the radius of the mirror is increased, the deviations re- 
tain their shape but decrease in size. For R = 20 cm, 
the beam is very close to the infinite mirror expectation. 
There is still a difference near the outer edge of the beam, 
causing the intensity to decay more quickly with radius 
than in the infinite mirror case. 

Our numerical results suggest that the anomalous 
diffraction loss is related to the deviation of the beam 
from the ideal Mesa shape. As D increases, the clipping 
approximation predicts a smooth increase in diffraction 
loss due to the widening beam. Finite mirror effects in- 
crease with the ratio D/R, so they alone do not explain 
this unexpected behavior. For the values of D that yield 
an anomalous diffraction loss, the variations around the 
plateau (Fig. ^}p) have an organized shape with an ap- 
proximate wavelength of wq. The variations in these 
cases have a shape such that they alter the falloff of the 
beam at the edge, i.e. in the last two centimeters of 
the mirror. This has an obvious beneficial effect on the 
diffraction loss. 

We expect that the fundamental mode of the cavity 
(the hyperboloidal shape, with no nodes), should have 
the lowest diffraction loss. However, this is not the case 
for all choices of a and D. Surprisingly, the diffraction 
loss of the first excited axisymmetric mode can decrease 
below that of the fundamental mode. This occurs, for 
a = TT, for a small range of D around 10.5 cm, with 
diffraction losses around 4 ppm. The arguments of two 
eigenvalues remain well separated. Fig.[5]shows the cross- 
ing of the two diffraction loss curves. 
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SO the diffraction loss of the fundamental mode is now 
lowest. Further increasing the mirror radius to 20 cm 
causes the losses of the second and third excited modes 
to cross. Despite changes of the mirror radius, the argu- 
ments of the eigenvalues change by less than a percent 
(we are only considering eigenvalues with losses less than 
10,000 ppm because otherwise they would dissipate too 
quickly to be of interest). This dependence on mirror 
radius tends to confirm that this is a finite mirror effect. 
Having demonstrated that this effect can have substan- 
tial and beneficial effects on the diffraction loss, we now 
turn our attention to studying the parameter space of 
hyperboloidal beams in more detail. Further studies of 
the precise cause of the anomalous losses may want to 
focus on the deviation of the beam from its theoretical 
expectation, as well as the behavior with changing mirror 
radius. 



B. Noises for fixed D 

The width of the beam increases with increasing D, 
which averages the fluctuations over more of the mirror 
surface and therefore decreases the noise. It is less clear 
how the noise will behave when a is changed. We be- 
gin by fixing D — cm and ignoring the diffraction 
loss constraint. Fig. [6] shows that all three types of noise 
increase as a moves away from and tt. Substrate ther- 
moelastic noise is most affected by changing a, followed 
by the coating noises (recall that both types of coating 
noise follow the same scaling law). 
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FIG. 5: Tfie diffraction losses of the fundamental mode and 
the first excited mode are plotted as a function of D for the 
a = TT Mesa configuration. The two curves cross due to finite 
mirror effects causing an anomalous diffraction loss for the 
first excited mode. 



FIG. 6: The noise ratios Noises /NoiseMcsa — 1 for three types 
of noise are shown as a function of a for fixed D = 10 cm. 
The minimal Gaussian a = 7r/2 has the highest noise and 
Mesa (q — 0, tt) has the lowest noise. 



To investigate the cause of this crossing, we choose the 
specific value D = 10.52 cm and increase the mirror ra- 
dius to 18 cm. The diffraction loss of the fundamental 
mode decreases from 3.7 ppm to 0.4 ppm, while the first 
excited mode only decreases from 3.6 ppm to 2.3 ppm 



As a is decreased from tt, the noises increase if D is 
kept fixed. At the same time, the diffraction loss de- 
creases. If we keep the diffraction loss fixed, the D can 
be increased as a decreases toward 7r/2. Widening the 
beam tends to decrease the noise, which partially offsets 
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the increase from changing a. We fixed the loss at 1.4 
ppm and found that the noise still increases for beams 
other than Mesa. We expect that for larger difi^raction 
losses this result will still hold. However, at 1 ppm the 
anomalous behavior of the diffraction loss due to finite 
mirror effects changes this conclusion. 



Hyperboloidal beams with 1 ppm Diffraction 
Loss 



As a decreases from tt toward n/2, the beam loses its 
flat top and sharp falloff, and approaches the minimal 
Gaussian. Also, the D corresponding to the local mini- 
mum in diffraction loss increases, and the local minimum 
becomes deeper and wider. Fig. [3] shows the diffraction 
loss versus D for three values of a, while Fig. [4] gives D 
and diffraction loss at the local minimum for a range of 
a. The local minimum of the diffraction loss for a Mesa 
beam {a = tt) is at D = 10.67 cm, and has 1.4 ppm. As 
discussed above, the Mesa is has the lowest noise of the 
hyperboloidal family for this diffraction loss. 

If a diffraction loss of strictly 1 ppm is required, the 
beam must be reduced to a width of D = 9.62 cm in the 
Mesa configuration. Alternatively, we can consider other 
values of a. The local minima are displayed in Fig. |4l 
which gives the values of D as well as the corresponding 
diffraction losses. Note that for a — 0.9tt, the diffraction 
loss at the local minimum is now below 1 ppm. Fig. [7^ 
shows the maximum D that yields a 1 ppm loss, as a is 
varied. There is a discontinuity between a = 0.9l7r and 
0.927r because below 0.927r the diffraction loss at the lo- 
cal minimum is below 1 ppm. The noises therefore drop 
substantially when a = 0.9l7r as in Fig.lTJj. The coating 
noise decreases by 12% and the substrate thermoelastic 
by 19%, relative to the 1 ppm Mesa beam. A strict re- 
quirement of 1 ppm diffraction loss therefore combines 
with the finite mirror effects to make the a = 0.9l7r, 
D — 10.94 cm configuration the best choice. 



D. Correcting for Finite Mirror Effects 

Restructuring the mirror to specifically account for 
finite mirror effects allows us to increase D in hyper- 
boloidal beams (thereby reducing thermal noise), while 
keeping within given diffraction loss constraints. The re- 
structured beams can reduce the diffraction loss by a fac- 
tor of 30 to 100, allowing a wider beam. For the Mesa 
case this allows for a net noise reduction of 30% for the 
beam satisfying the Ippm diffraction loss constraint. 

As noted above, the original Mesa beam used to con- 
struct the mirror is infinite in extent. The mirror is de- 
signed to be a phasefront of the theoretical beam. Since 
the mirror is actually finite, for D ^ R there can be sub- 
stantial effects due to missing light that was incident on 
the mirror plane outside the mirror radius. To account 
for these finite mirror effects we reconstruct the mirror, 
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FIG. 7: (a) The largest values of D giving a strict 1 ppm 
diffraction loss are shown as a function of a. The discontinuity 
is caused by finite mirror effects as explained in the text, (b) 
The corresponding values of the noise, normalized so that the 
noises for the 1 ppm Mesa {a = tt, D — 9.62 cm) are all equal 
to 1. The discontinuity in allowed maximum D leads to a 
sharp drop in noise at q = O.QItt. 



with the goal of making the phase of the first eigenbeam 
constant at the mirror surface, rather than the phase of 
the idealized infinite beam. 

The propagator formulation allows us to explicitly cal- 
culate the phase of the eigenbeams, as a function of r. As 
the mirror deviation from the fiducial spheroid h{r) en- 
ters into the calculation of phase through the propagator, 
we use an iteration scheme to adjust the mirror to match 
the eigenbeam phasefront motivated by the argument of 
the propagator: 



arg[/C(r, r')] 



7r/2 + fco ih{r) + h{r') - L) (26) 



_fc„cosa(- + -j 



We see in Eq. (pS)) that if the phase of the eigenbeam 
is too large at some radius r, reducing the value of h{r) 
should act to roughly reduce the phase of the new eigen- 
beam. 

With this motivation, we apply the simple iteration 
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FIG. 8: The phase of the fundamental eigenbeam as a function 
of radius is shown for iteration zero and 250. It can be seen 
that the iteration scheme drives it closer to zero as expected. 



(27) 



scheme: 

/i,+i(r) = h^ir) - c X arg [v^^^\r)/v}^\o) 

where the i denotes the i^^ iteration, V^'"^\r) is the first 
eigenbeam for the mirror with deviation hi{r), and c > 
is an arbitrary constant less than unity, chosen to prevent 
overshoot. 
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FIG. 9: The diffraction loss is shown as function of itera- 
tion number for a a = vr Mesa configuration with D = 11.35 
cm. The original beam (not shown) has diffraction loss at 
46.5 ppm, and it can be seen that after a few iterations 
the diffraction loss begins to converge to an exponential with 
lower bound ~ 1 ppm. The best fit exponential is given by 
0.96 + 1.616 exp(—0.013i) ppm, where i is the iteration num- 
ber. 



This iteration scheme successfully reduces the relative 
phase of the eigenbeam, as shown in Fig [51 Optimizing 
the mirror surface to match the phasefront of the pri- 
mary eigenbeam also acts to reduce the diffraction loss 
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Dilfraction Loss 
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a. = TV 


D = 11.35 


cm 
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Diffraction Loss 


Phase 




1.0 ppm 


0.0 




320.6 ppm 
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1100.1 ppm 
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a = O.gvr 


D = 11.87 


cm 




(Uniterated) 




Dilfraction Loss 


Phase 




43.5 ppm 


0.0 




205.0 ppm 


-0.4173 




371.2 ppm 


-1.1665 






-2.1082 




"1072^ fi nnm 


3.0878 




a = 0.97r 


D = 11.87 


cm 




(Iterated) 




Dilfraction Loss 


Phase 




1.0 ppm 


0.0 




366.6 ppm 


-0.4185 




2013.4 ppm 


-1.1701 




32628.3 ppm 


-2.1267 




87359.8 ppm 


3.0335 





TABLE III: The phase separation for the axisymmetric (m = 
0) modes with diffraction losses less than 10%, both before 
and after the iteration scheme is applied. The phases do not 
change significantly as the mirror is iterated. The absolute 
value of the eigenvalues with nonzero phase (and hence the 
diffraction loss) increase upon iterating. Preliminary calcu- 
lations show that non-axisymmetric modes have diffraction 
losses increased by the iteration process while the phases 
change by no more than 5%. 



for that mode in general, with the iteration scheme pro- 
viding convergence towards an apparent lower bound for 
the diffraction loss, while increasing the diffraction loss 
for other higher-order eigenbeams. This lower bound in- 
creases with D fFig [TO)l . 

This diffraction loss is plotted against the iteration 
number for the Mesa {a = tt) case with D — 11.35 cm in 
FigUl The iteration scheme is shown to lower the diffrac- 
tion loss for this D from 46.5 ppm to a ~ 1 ppm lower 
bound, satisfying the required design constraint. The 
diffraction losses of higher order modes are more than 
doubled in the iterated case as illustrated in Table [1111 

The beam for the iterated mirror with D — 11.35 cm 
is close to the original Mesa, but with variations in the 
central plateau of relative amplitude ^ 1/30 and varia- 
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a 


D [cm] 


Coating 


Substrate 


Substrate 






Noise 


Brownian 


Thermoelastic 








Noise 


Noise 




11.35 


0.72 


0.84 


0.63 


0.9n 


11.87 


0.80 


0.90 


0.69 



TABLE IV: The coating, substrate Brownian and substrate 
thermoelastic noise are displayed after the iteration process. 
The diffraction loss is kept constant at 1 ppm. The noises 
are normalized to noises of the original (a = n) Mesa (with 
D — 9.62 cm) that gives the 1 ppm diffraction loss. It can be 
seen that the iteration scheme lowers the noise by about 30% 
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FIG. 10: The diffraction loss in ppm calculated using the clip- 
ping approximation is compared to that using the propagator 
eigenvalues for the iterated and original mirrors as a function 
of D. As before, the configuration studied is a = tt Mesa with 
i? = 17 cm. The iteration process lowers the diffraction loss 
by a factor of 30 to 100. 



tions of radial wavelength ~ wq fFig fTTjl . This seems to 
be an unavoidable consequence of a finite R, as even the 
original eigenbeam has roughly similar features. Despite 
changing the variations in the plateau of the beam inten- 
sity the process of iteration does not significantly affect 
the noises computed using Eq. 

Similarly the iterated mirror has variations of the cen- 
tral mirror shape of similar radial scale, with amplitude 
on the order of 2 nm, shown in Fig [121 The most signif- 
icant feature of the iterated mirror is the inward tilting 
of the outer edge of the mirror, though preliminary stud- 
ies show that the primary contribution to reducing the 
diffraction loss is due to the variations near the center of 
the mirror. 

Reformulating the mirror to account for finite mirror 
effects allows us to increase the D parameter of the beam 
from 9.62 cm to 11.35 cm for a concentric Mesa beam 
while still maintaining a Ippm diffraction loss. Though 
this design may introduce more complications in the con- 
struction of the mirror itself, it allows a significant re- 




r [cm] 



FIG. 11: The intensity profile \Utt\ for the mirror with R = 
oo, i? = 17 cm uniterated and R — 17 cm at iteration 250 
are compared for the a = n Mesa configuration with D = 
11.35 cm. The finite mirror effects induce oscillations in the 
intensity profile that do not disappear when the mirror is 
corrected. An inset shows the central 8 cm 'plateau' of the 
beam in detail. 
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FIG. 12: The correction to the mirror Htt at iteration and 
at iteration 250 are compared for the a = tt Mesa configura- 
tion with D — 11.35 cm, R — 17 cm. The iteration scheme 
introduces some bumps on the mirror of the size ~ 2 nm. The 
inset shows the central 8 cm of the mirror in more detail. 



duction in noise by broadening the beam. This iteration 
scheme can also be used for other values of a, as shown 
in Table III, where the iterated mirror for a = O.Ott has 
a diffraction loss lower bound of 1 ppm for D = 11.87 
cm. However, we find that a = tt is optimal for noise 
reduction. (See Table [IVl) 



VII. CONCLUSIONS 

In this paper, we studied thermal noise and diffrac- 
tion loss for the hyperboloidal family of light beams and 
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mirror shapes in detail for the first time. This fam- 
ily had been initially proposed to unify the concentric 
{a — tt) and nearly-flat {a = 0) Mesa configurations 
through variations of the twist angle a. In this paper 
we also presented an analytic proof that Gaussian beams 
are a limiting case of the hyperboloidal beam as I? — > oo. 
This was previously conjectured in Ref. [l3|- We de- 
veloped a pseudo-spectral code both fast and accurate 
enough to calculate the diffraction loss directly from the 
beam propagator. We find that the finite radius of the 
mirror causes beam shapes to deviate significantly from 
the infinite-mirror theoretical expectations. This causes 
a previously unnoticed local minima in the diffraction 
loss that can be exploited to find a natural beam width 
D for the current diffraction loss constraints of about 1 
ppm. For an a = TT Mesa beam a local minimum oc- 
curs at D = 10.67 cm giving 1.4 ppm diffraction loss for 
a mirror of radius R = 17 cm. If one requires a strict 
enforcement of the 1 ppm diffraction loss we show that 
a hyperboloidal beam with a — 0.9l7r and D ~ 10.94 
cm has lower noise than that of the a = tt Mesa with 1 
ppm diffraction loss. The coating noise changes by about 
12% and the substrate thermoelastic noise and substrate 
Brownian noise change by 5% and 19%, respectively. 

We also propose new mirror and beam shape config- 
urations that explicitly account for finite mirror effects 
by reformulating the mirror surface to coincide with the 
phasefront of the primary eigenbeam. These beams re- 
duce the diffraction loss by more than an order of mag- 
nitude for the range of D considered here (between 10 
cm and 11.8 cm). This allows the use of wider beams 
for the same diffraction loss constraints on the primary 
eigenmode, while the diffraction losses of higher order 
modes (both axisymmetric and non-axisymmetric) are 
increased. We are able to widen the a = tt Mesa beam 
for a i? = 17 cm mirror while keeping the diffraction loss 
fixed at 1 ppm from a width oi D = 9.62 cm to Z? = 11.35 
cm. This lowers the coating thermal noise by about 30% 
(compared to the smaller D Mesa) and the other noises 
(substrate Brownian and substrate thermoelastic noise) 
by comparable factors. However, feasibility of the con- 
struction of the mirrors must also be taken into account. 

The non-iterated beams taking advantage of the lo- 
cal minimum in diffraction loss discussed above are sup- 



ported by strict hyperboloidal or Mesa mirrors, which 
would be no harder to make than the current Mesa de- 
signs, and would still lower the coating thermal noise by 
12%. If one is to consider the more ambitious goal of 
lowering the coating thermal noise by 28% while using a 
beam that is very similar to Mesa through the iteration 
scheme described, the limitations on mirror manufactur- 
ing errors are likely to be more stringent, but still less 
than the currently considered conical beams [1, [l9[ . In 
addition, the methods developed here for reducing the 
diffraction loss of the Mesa beam may be applied in the 
case of the conical beams previously considered [1, [l^ . 
The phase fronts of conical beams considered there have 
not been optimized to match the finite mirror surface. 

Recently, parametric instability [2l|, [23, H^, was 
found to be a serious problem in Advanced LIGO. Choos- 
ing D at the minimum of the diffraction loss curve of the 
hyperboloidal beams increases the diffraction loss of the 
higher eigenmodes, thus in principle somewhat improv- 
ing the parametric instability. The effect is most pro- 
nounced in the case of the iterated mirrors which also 
have the most drastic thermal noise reduction. 
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